EQUALIZER FOR COMMUNICATION OVER NOISY CHANNELS 



FIELD OF THE INVENTION 



The present invention relates generally to data communication systems and more 



particularly relates to an apparatus for and a method of iterative equalization for 
5 communication over noisy channels. 



In recent years, the world has witnessed explosive growth in the demand for all types 
of communications and it is predicted that this demand will increase in the future. In order to 
communicate information accurately and efficiently, communications receivers are designed 
10 to extract information from the channel as best they can. Receivers are designed to operate 
optimally with certain types of channels, e.g., power line, IR, RF, optical, twisted pair, coax, 
3 etc. Various types of noise, however, are a characteristic feature of most channels. Well- 
;n known types of channel noise include Intersymbol Interference (ISI) caused by pulse 
,4 spreading in the channel, fading, intermodulation distortion, multipath, etc. Performance is 
:: 15 severely affected in receivers not designed to adequately handle such types of channel noise. 

Equalization is a well-known signal processing technique used to combat the 
!;=;: Intersymbol Interference distortion imparted to the transmitted by the channel whereby the 
4 ; receiver attempts to compensate for the effects of the channel on the transmitted symbols. An 
j equalizer attempts to determine the transmitted data from the received distorted symbols 
using an estimate of the channel that caused the distortions. Examples of commonly used 
types of equalizers include the maximum likelihood sequence estimation (MLSE) equalizer 
that utilizes the well known Viterbi Algorithm (VA), linear equalizer and decision feedback 
equalizer (DFE). In communications systems where ISI arises due to partial response 
modulation or a frequency selective channel, a maximum likelihood sequence estimation 
25 (MLSE) equalizer is optimal. 

Many modern equalization techniques, however, utilize very complex signal 
processing algorithms to achieve acceptable levels of performance. Practically, most of these 
techniques can only be performed on expensive high powered digital signal processors. With 
the ever pressing demand to make communications enabled products smaller, cheaper and 
30 high- performing, it would be desirable to implement the entire communications hardware in 
silicon, such as in an Application Specific Integrate Circuit (ASIC) or the like. To place the 
complex, sophisticated equalization techniques of the pri^^^a|^^^^ 5j ^Qwever, would 
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be nearly impossible in terms of gate count and cost, given today's available gate densities 
and chip sizes. Prior art equalization techniques require extremely large processing resources 
and memory in order to implement them in silicon. 

Thus, there is a need for an equalization technique that is capable of achieving 
adequate levels of performance while lending itself for efficient and inexpensive 
implementation in silicon, e.g., ASICs, etc. 
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SUMMARY OF THE INVENTION 
The present invention is an apparatus for and a method of iterative equalization for 
communication over noisy channels. The present invention is suitable for use with a wide 
range of different types of communications systems and channels and is particularly 
5 applicable to power line based communication systems. 

Many types of channels, especially the power line media, typically suffer from 
Intersymbol Interference due to the pulse spreading effects of the channel. The iterative 
equalization technique of the present invention provides a solution to the communication 
problems imposed by the characteristics of the power line media, e.g., Intersymbol 
10 Interference, etc. In addition, the technique is well suited for implementation in integrated 
circuit form, e.g., ASIC, etc. thus enabling the communications over noisy channels in 
receivers employing a wide range of communication technologies. 
3 The detection of the received bits is performed by first eliminating the influence of the 

; : : channel hi using the iterative equalization technique of the present invention. This technique 
115 comprises two equalization steps: the first step performs iterative subtraction while the 
second step uses the inverse of the transfer matrix. Note that the iterative equalization 
- ;| method of the present invention may be applied to any digital modulation, such as BPSK, 
QPSK and QAM modulations. Although the invention is described in the context of BPSK 
symbols, it is appreciated that one skilled in the art can apply it to communication systems 
20 wherein each symbol represents a plurality of bits. 

u In iterative subtraction, the impulse response of the best estimate of previously 

estimated symbols (which include both past and future symbols) are processed by a soft 
decision function. The results are then subtracted out to yield the impulse response of a 
plurality of adjacent received symbols, termed the modified received symbol. This operation 

25 'cleans' the intersymbol interference due to the previously processed symbols with the 
exception of a few neighboring symbols. Once the modified received symbol is calculated, 

the estimate of the current symbol B n can then be calculated by passing the modified 

received symbol through a matched filter which functions to collect the energy within the 
symbol. 

30 Note that throughout this document, the following definitions apply. The term 

'previous symbol' is defined as a symbol that has already been processed and whose index is 
lower than the current symbol being processed by the equalizer, i.e. B n where n>0. The index 
refers throughout to the time when the symbol was sampled by the A/D converter. The term 
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'future symbol' is defined as a symbol that is received later in time and that has an index 
higher than that of the current symbol. The term 'previously processed symbol' is defined as 
a symbol that is still within the equalizer and may comprise a previous or future symbol with 
regard to its index. 

For example, consider the iterative equalizer in iteration j and processing symbol n. 
In this example, the previous symbols, indexed n-l...n-m, were last processed during iteration 
j. Future symbols, indexed n+l...n+m were last processed during iteration y-L Note that the 
term 'previously estimated symbol' is identical to the term 'previously processed symbol.' 

It is important to note that bad estimates of the previously processed symbols may 
influence the estimation of the current symbol. In order to overcome this problem, the 
present invention performs one or more iterations. In the second iteration, the majority of the 
influence from neighboring symbols is already subtracted out, thus enabling more accurate 
estimations to be calculated. Third and subsequent iterations are performed to achieve even 
more accurate results. 

In the example embodiment presented herein, the influence of the n-\ th previously 
processed symbol is calculated using the inverse matrix derived using the optimum 
equalization technique rather than using iterative subtraction. This enhances the accuracy of 
the resultant soft symbol decision. 

There is therefore provided in accordance with the present invention a method of 
equalization of symbols transmitted over a channel having a channel impulse response hi, the 
method comprising the steps of calculating the impulse response of previously estimated 
symbols utilizing the channel impulse response hi and previously estimated decisions, 
selecting a window having a size of one or more symbols and removing the influence of 
previously estimated symbols from the received signal so as to generate a modified received 
signal, correlating the modified received signal with the channel impulse response to obtain a 
correlation for each symbol in the window, estimating the results of the correlation so as to 
yield a new estimate of one or more of the transmitted symbols and repeating the steps of 
calculating, selecting, removing, correlating, estimating and updating. 

There is also provided in accordance with the present invention a method of 
equalization of symbols transmitted over a channel having a length m and a channel impulse 
response hi, the method comprising the steps of providing a sample shift register array 
comprising N rows wherein N represents the number of samples per symbol, calculating a 
cross correlation coefficient matrix and an inverse matrix determinant based on a channel 
estimate, generating the impulse response of a first number of previous symbols, subtracting 
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the impulse response of the first number of previous symbols from the contents of the register 
array, cyclically shifting the contents of the register array through a second number of cycles 
constituting an iteration and repeating the steps of generating and subtracting so as to 
generate a modified received symbol wherein the influence of the first number of previously 
5 processed symbols is removed, performing a matched filter operation on the contents of the 
shift register array, generating an estimate of a symbol in accordance with the output of the 
matched filter, removing the influence of previously estimated symbols using the previously 
estimated symbols and the cross correlation coefficient matrix so as to yield symbol output B n 
and wherein m and N are positive integers. 
10 There is further provided in accordance with the present invention an equalizer for 

equalizing, over a plurality of iterations, symbols transmitted over a channel having a channel 
impulse response hi comprising first means for iteratively removing the influence of previous 
; I estimated symbols from a received symbol so as to generate a modified current received 
: :: symbol and for correlating the modified current received symbol with the channel response to 
: X5 yield a matched filter output, a new modified received symbol and matched filter output 
rj generated every iteration of the equalizer, second means for calculating a current symbol 
decision utilizing a previous symbol decision, the output of the matched filter and a cross 
correlation coefficient matrix and a quantizer adapted to quantize the current symbol decision 
"r to a quantized symbol decision having a finite number of possible values and third means for 
&) calculating a correction value comprising the difference between the current symbol decision 
; ~: calculated during the current iteration and the previous symbol decision calculated during the 
previous iteration. 

There is also provided in accordance with the present invention a communications 
receiver for receiving and decoding an M-ary transmitted signal transmitted over a channel 

25 having a channel impulse response ly comprising a front end circuit for receiving and 
converting the M-ary transmitted signal to a baseband signal consisting of a plurality of 
received samples, a demodulator adapted to receive the received samples and to generate 
input samples therefrom in accordance with the M-ary modulation scheme used to generate 
the transmitted signal, an iterative equalizer operative to receive the input samples and to 

30 generate a sequence of symbol decisions therefrom, the iterative equalizer comprising 
processing means programmed to: calculate the impulse response of previously estimated 
symbols utilizing the channel impulse response hi and previously estimated decisions, select a 
window having a size of one or more symbols and remove the influence of previously 
estimated symbols from the received signal so as to generate a modified received signal, 
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correlate the modified received signal with the channel impulse response to obtain a 
correlation for each symbol in the window, estimate the results of the correlation so as to 
yield a new estimate of one or more of the transmitted symbols, repeat the steps of 
calculating, selecting, removing, correlating, estimating and updating, a decoder adapted to 
5 receive the estimate and to generate binary received data therefrom and wherein M is a 
positive integer. 

There is also provided in accordance with the present invention a computer readable 
storage medium having a computer program embodied thereon for causing a suitably 
programmed system to equalize symbols transmitted over a channel having a channel impulse 
10 response hi by performing the following steps when such program is executed on the system: 
calculating the impulse response of previously estimated symbols utilizing the channel 
impulse response hi and previously estimated decisions, selecting a window having a size of 
i one or more symbols and removing the influence of previous symbols from the received 
:V; signal so as to generate a modified received signal, correlating the modified received signal 
15 with the channel impulse response to obtain a correlation for each symbol in the window, 
i estimating the results of the correlation so as to yield a new estimate of one or more of the 
\[: transmitted symbols, repeating the steps of calculating, selecting, removing, correlating, 
estimating and updating and wherein M is a positive integer. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
The invention is herein described, by way of example only, with reference to the 
accompanying drawings, wherein: 

Fig. 1 is a block diagram illustrating an example communications system including 
5 transmitter, channel and concatenated receiver; 

Fig. 2 is a block diagram illustrating the structure of an example concatenated receiver 
including an iterative equalizer constructed in accordance with the present invention; 

Fig. 3 is a flow diagram illustrating the iterative equalization method of the present 
invention; 

10 Fig. 4 is a block diagram illustrating the iterative equalizer of the present invention in 

more detail; 

Fig. 5 is a block diagram illustrating the modified received bit calculator/matched 
filter of Figure 4 in more detail; and 
\ :l Fig. 6 is a block diagram illustrating an example computer system adapted to perform 

U-5 the iterative equalization method of the present invention. 
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DETAILED DESCRIPTION OF THE INVENTION 
Notation Used Throughout 



The following notation is used throughout this document. 



Term 


Definition 


ASIC 


Application Specific Integrated Circuit 


BER 


Bit Error Rate 


BPSK 


Binary Phase Shift Keying 


CPU 


Central Processing Unit 


DEE 


Decision Feedback Equalizer 


DSP 


Digital Signal Processing 


EEROM 


Electrically Erasable Read Only Memory 


EEC 


Forward Error Correction 


FPGA 


Field Programmable Gate Array 


IEEE 


Institute of Electrical and Electronic Engineers 


IR 


Infrared 


ISI 


Intersymbol Interference 


LMS 


Least Mean Squares 


LSE 


Least Squared Error 


MAP 


Maximum A Priori 


MLSE 


Maximum Likelihood Sequence Estimation 


MMSE 


Minimum Mean Squared Error 


PN 


Pseudo Random 


QAM 


Quadrature Amplitude Modulation 


QPSK 


Quadrature Phase Shift Keying 


RAM 


Random Access Memory 


RF 


Radio Frequency 


ROM 


Read Only Memory 


SOVA 


Soft Output Viterbi Algorithm 


VA 


Viterbi Algorithm 



Detailed Description of the Invention 

The present invention is an apparatus for and a method of iterative equalization for 
communication over noisy channels. The present invention is suitable for use with a wide 
range of different types of communications systems and channels and is particularly 
applicable to power line based communication systems. 

Many types of channels, especially the power line media, typically suffer from 
Intersymbol Interference due to the pulse spreading effects of the channel. The iterative 
equalization technique of the present invention provides a solution to the communication 
problems imposed by the characteristics of the power line media, e.g., Intersymbol 
Interference, etc. In addition, the technique is well suited for implementation in integrated 
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circuit form, e.g., ASIC, etc. thus enabling the communications over noisy channels in 
receivers employing a wide range of communication technologies. 

A block diagram illustrating an example communications system including 
transmitter, channel and concatenated receiver is shown in Figure 1. The communications 
system, generally referenced 100, employs an inner and outer encoder in the transmitter and 
inner and outer decoding stages in the receiver. The communication system 100 represents 
the typical scheme that may be used in systems that utilize a wide range of different channel 
media, e.g., power line, IR, RF, optical, coax, twisted pair, etc. In such a system, the 
transmitter 102 comprises an encoder 106, interleaver 108, bit to symbol mapper or symbol 
generator 110, modulator 112 and power amplifier/coupling circuitry 114. Input data bits 
104 to be transmitted are input to the encoder 106, which may comprise an error correction 
encoder such as a Reed Solomon encoder, a convolutional encoder, a parity bit generator, etc. 
The encoder functions to add redundancy bits to enable errors in transmission to be located 
and fixed. 

In the implementation of the iterative equalizer presented herein, two assumptions are 
made: (1) the channel parameters are constant for the duration of the entire packet, and (2) 
the channel response for the transmitted symbol stream is the sum of the channel response for 
each symbol, or in other words, that the channel is linear. In practice, channel parameters can 
be tracked during the packet using known adaptive schemes. 

It is noted that the channel 120 often times can be modeled as an L-symbol long FIR- 
type channel, and the encoder and the channel can be viewed as outer and inner encoders. 
Both the inner and outer encoders in the transmitter have complementary decoders in the 
receiver. 

The bits output of the encoder are then interleaved by interleaver 108. Interleaving 
rearranges the order of the bits so as to help combat burst errors encountered in the channel or 
burst errors resulting from the operation of the equalizer. The bits output of the interleaver 
108 are then mapped to symbols by the bit to symbol mapper 1 10. The bit to symbol mapper 
functions to transform the bits to modulator symbols. 

The output from the mapper 110 is input to the modulator 112 which receives 
symbols in the M-ary alphabet and generates the analog signal that is subsequently 
transmitted over the channel 120 via the transmit amplifier/coupling circuitry 1 14. Note that 
the channel may comprise a wired channel, e.g., power line, twisted pair, xDSL, ISDN, 
Ethernet, etc. or may comprise a mobile wireless channel, e.g., cellular, cordless, satellite, 
etc. In each case, the processing performed in the transmitter is intended to generate a signal 
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that can be transmitted over the channel so as to provide robust, error free detection by the 
receiver. 

At the receiver 122, the analog signal from the channel is input to front end circuitry 
124 which demodulates and samples the received signal to generate received samples x(k) 
132. The symbols are first input to an inner decoder 126. An example of an inner decoder is 
an equalizer which compensates for the ISI caused by the delay and time spreading of the 
channel. The function of the equalizer is to detect the symbols that were originally 
transmitted by the modulator. Examples of commonly used types of equalizers include the 
maximum likelihood sequence estimation (MLSE) equalizer that utilize the well known 
Viterbi Algorithm (VA), linear equalizer and decision feedback equalizer (DFE). Examples 
of soft output type equalizers include Soft Output Viterbi Equalizers (SOVA) and Maximum 
A Posterior (MAP) algorithm based equalizers. 

The output of the inner decoder comprises symbols s(k) 134 that represent either hard 
or soft decisions. The decisions are then input to a de-interleaver 128. The de-interleaved 
stream is then input to an outer decoder 130 which functions to locate and fix errors using the 
redundancy inserted by the encoder. The outer decoder generates the binary receive data a(k) 
138. 

An example of an outer decoder is a convolutional decoder that utilizes the Viterbi 
Algorithm. The Viterbi algorithm is widely used in communication systems and has been 
adapted to perform functions including demodulation, decoding, equalization, etc. Many 
systems utilize the Viterbi Algorithm in both the inner and outer decoding stages. Other 
examples include diagonal parity and even odd parity based decoders. 

The outer decoder, in some systems, is adapted to utilize the hard decisions output 
from the inner decoder, e.g., the equalizer. Optimal decoders, on the other hand, require soft 
decisions instead. For example, an outer decoder that utilizes the Viterbi Algorithm to 
perform convolutional forward error correction decoding, requires soft decisions as input. 
The advantage of a Viterbi decoder is that it can efficiently process soft decision information. 

Equalization is a well known technique used to combat intersymbol interference 
whereby the receiver attempts to compensate for the effects of the channel on the transmitted 
symbols. An equalizer attempts to determine the transmitted data from the received distorted 
symbols using an estimate of the channel that caused the distortions. Stated in other words, 
equalization is a technique use for removing deterioration caused by the channel from 
received symbols, wherein the deterioration is usually caused by intersymbol interference. 
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The ISI may be caused by reflections, multipath, fading and the effects of a non-flat transfer 

function for the channel. 

Mathematically, the problem can be described as transmitting a series of symbols (or 

bits) through a channel having an impulse response denoted by Note that throughout this 
5 document, the transmission of either symbols or bits holds regarding application of the 

present invention. The invention applies to the general case of M-ary signaling, the binary 

case being when M=2. We will continue with the binary case as an illustrative example in the 

following description. It is appreciated that the description below can be applied to the 

general case of M-ary signaling. 
10 For simplicity and clarity sake, it is assumed that the receiver input is sampled at a 

rate of one sample per symbol. All the equations presented infra can be modified by one 

skilled in the art to reflect the case of fractional spaced equalization wherein the receiver 
~Z input is sampled at a rate of more than one sample per symbol. The output of the channel x n 

can be expressed as 

*5 x^Wn+J^K-i-hi (1) 

V[ where 

x n represents the received samples 
7: B n represents the bits transmitted (e.g., +1, -1); 

w n represents the noise; 
20 h n represents the impulse response of the channel; 

m represents the length of the channel, i.e. the length of the impulse response 

In accordance with the invention, the detection of the received bits is performed by 
first eliminating the influence of hi using the technique of iterative equalization described 
25 herein. This technique comprises two equalization steps: the first step performs iterative 
subtraction while the second step uses the inverse of the transfer matrix. Both methods are 
described in more detail hereinbelow. 

First Iteration 

The method assumes that the impulse response hi is either known or estimated by 
30 some external means, e.g., correlation with a known preamble with good autocorrelation 
properties, and that the estimates of the previous symbols (bits) B n _ x ,..B n _ m were previously 
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determined, wherein m is the length of the impulse response. Then the following expression 
can be used to calculate the impulse response of the current symbol (bit) 

m 

yj =x j- H B h-K J=n...n+m-l (2) 

where 

yj represents the modified received samples with the influence of all previous symbols 
removed; 

xj represents the received symbol (or bit); 

B n represents the bits transmitted; 

h n represents the impulse response of the channel; 

m represents the length of the channel, i.e. the length of the impulse response 

Note that the entity yj represents the impulse response of the symbol n (including additive 
noise) only whereby the impulse response of the previous symbols has been subtracted out. 
This operation 'cleans' the intersymbol interference from previous symbols but not from 
future symbols, so still a large amount of ISI is left. Note that, for purposes of calculating 
Equation 2, the soft decisions (as are computed below) are used for the previous symbols 
while zero is used for the case j-i<0 

Once the modified received symbol y } is calculated, the estimate of the current symbol 
(bit) B n can then be calculated using the following expression 

m 

5„=/(£j i+ „-V) (3) 

i=0 

where j{ ) is a soft decision function (e.g., hyperbolic tangent). The slope of the decision 
function preferably increases with the iterations. Thus, the modified received symbol y„ is 
processed through a matched filter operation (i.e. a correlation function) which functions to 
estimate the symbol. Note that in the case of QPSK the soft decision is applied separately on 
the real and imaginary component. In the case of QAM or other modulation, the soft decision 
needs to be modified accordingly. 

It is important to note that bad estimates of the previous symbols (bits), and the ISI 
due to future symbols may influence the estimation of the current bit. In order to overcome 
this problem, the present invention performs one or more iterations. In addition, the matrix 
inversion method as described below is used to increase the performance of the equalizer. 
After the second iteration, the majority of the influence from neighboring symbols has 
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already been subtracted out, thus enabling more accurate estimations to be calculated. Third 
and subsequent iterations are performed to achieve even more accurate results. During the 
additional iterations, the future symbols (i.e. symbols indexed n+1 .. n+m) are estimated and 
subtracted as is shown infra. 

The iteration process can be performed a fixed number of times or can be repeated 
until the Bit Error Rate (BER) or some other metric converges to a constant (i.e. a minimum 
error is achieved in the estimate). Note that since the influence of the neighboring symbols 
was already subtracted from the receive signal, only the error of the neighboring symbols 
need be calculated and compensated in subsequent iterations (after the first iteration). 

Improvement by Transfer Matrix Calculation 

In accordance with the iterative equalization method of the present invention, the 
output of the matched filter, can be expressed using the following expression 

m 

*»=E*rf«'V (4) 

where 

z n represents the output of the matched filter; 

x n represents the input to the equalizer; 

h n represents the impulse response of the channel; 

m represents the length of the channel, i.e. the length of the impulse response 

Note that the format of Equation 4 is similar to that of Equation 3. Equation 4 
implements a matched filter. The relationship between the transmitted bits B n and the output 
of the matched filter z n neglecting the noise can be expressed as shown in Equation 5. 

Note also that for the complex case, where both h and x{n) are complex, the upper 
triangle in the matrix should be conjugated. It is appreciated that one skilled in the art would 
know how to modify the expression below for the complex case. 
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z n represents the output of the matched filter; 

matrix h is the channel estimate; 

B n represents the symbols (bits) transmitted; 

and 



i=0 

m 

^o=S fc i'Vi (6) 

m 



In conventional linear equalization, the approach to determine B n is to calculate 
matrix h and then compute its inverse. This approach, however, is very difficult to 
implement in an ASIC and is not practical. This is because the calculation of the inverse of 
the matrix is very difficult when implementing this algorithm in hardware or software. 

Thus, in accordance with the present invention, a relatively small number of 
neighbors are calculated using the inverse matrix equalizer. All the other neighbors are 
cancelled out using the iterative subtraction equalization technique of the present invention. 
Note that any number of neighbors may be calculated depending on the processing resources 
available. Through experimentation, it has been found that calculating just one or two 
neighbors yields good performance. Note that any linear equalization could be used 
including, but not limited to, zero forcing, Least Mean Squares (LMS), matrix inversion 
using the MMSE criterion and matrix inversion using the LSE criterion. 

In the next iteration, Equation 2 is modified to yield 

m j -n- S 

yj= x j- Z S h h i-T B H ■ h, , j=n. . .n+m-1 (7) 

i=j~n+\ t=Q 

where the third term in the equation represents the contribution of future symbols to x. j and 

S represents the number of neighbors to be processing using linear equalization. These S 
symbols are making up a sliding window. 

The output of the matched filter as applied to the y jf output of (7), neglecting the noise, 
can be expressed using the following expression 
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Z n yi+n * h 

Assuming the subtracted symbols are correct, we obtain 
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where the soft decision function is also applied to each element. If the window of S symbols 
is now slid one symbol to the right, only the last symbol from the estimated vector B n is 
needed for the feedback, so only this symbol needs to be actually calculated. In alternative 
embodiments of the invention, the window is slid by a number of symbols each time, and a 
number of symbols are calculated during the linear equalization step. 

A block diagram illustrating the structure of an example concatenated receiver 
including an iterative equalizer constructed in accordance with the present invention is shown 
in Figure 2. The equalization technique of the present invention is described in the context of 
a communications receiver. Note the communication receiver shown herein is presented for 
illustration purposes only and is not meant to limit the scope of the present invention. It is 
appreciated that one skilled in the communication arts can construct numerous other types of 
receivers incorporating the iterative equalization technique of the present invention. 

The receiver, generally referenced 10, comprises coupling circuitry 12 connected to 
the channel media, for example the power line, an A/D converter 14, synchronization/channel 
estimation 16, iterative equalizer 18, matrix calculation module 26, outer decoder 20, de- 
interleaver 19. 

The signal from the channel is input to the channel coupling circuit 12 which 
functions as the analog front end for the receiver. The channel coupling circuit comprises an 
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anti-aliasing filter (not shown). The signal output of the channel coupling circuitry 12 is 
input to the A/D converter 14. The digitized receive data stream is input to both the iterative 
equalizer 18 and the synchronization/channel estimation 16. Additional filtering, 
conditioning, decimation or down-conversion of the signal can optionally be applied before 
the equalizer. This module functions to detect the training sequence or sync sequence sent by 
the transmitter. This module also functions to generate an estimate of the channel. 

The channel estimation 16 is operative to generate a channel estimate 28 represented 
by hi that is used by both the matrix calculator 26 and the iterative equalizer 18. The channel 
estimation is generated using the received input samples x(n) 25. 

Numerous methods of channel estimation are known in the art and suitable for use 
with the present invention, for example, those based on correlation techniques and those 
based on least squares techniques. 

The output of the iterative equalizer 18 comprises soft symbol decisions B n 22. The 
iterative equalizer is operative to output soft decision information given (1) received signal 
from the channel x(n) 25, (2) channel estimate information h t 28 (e.g., FIR filter taps used to 
estimate the channel) and (3) channel estimate matrix element h w 27 The soft decision 
information B n output of the iterative equalizer 18 is input to the outer decoder 20 which is 
preferably an optimal soft decoder. The outer decoder functions to locate and fix errors using 
the redundancy bits inserted by the encoder. The outer decoder generates the binary receive 
data. Examples of the outer decoder 20 include convolutional decoders utilizing the Viterbi 
Algorithm, convolutional Forward Error Correction (FEC) decoders, turbo decoders, etc. 
Soft input Viterbi decoders have the advantage of efficiently processing soft decision 
information and providing optimum performance in the sense of minimum sequence error 
probability. 

Note that optionally, an interleaver/de-interleaver may be added to the system. In this 
case, a symbol based interleaver/de-interleaver is used. If a bit based interleaver/de- 
interleaver is used, some mechanism of mapping soft symbols to bits must be used before the 
outer decoder. 

In operation, the transmitter transmits a modulated signal to the channel. In the 
example embodiment described herein, the transmitter transmits a binary phase shift keying 
(BPSK) or quadrature phase shift keying (QPSK) modulated signal to the channel. The 
iterative equalization technique requires knowledge of the impulse response of the channel h t . 
This can be obtained using the channel estimation techniques described above. One method 



12407.0014 



16 



of determining hi is to transmit a training sequence or sync pattern comprising a pseudo noise 
(PN) sequence. The signals are transmitted grouped together in the form of packets wherein 
each packet begins with the PN synchronization sequence. 

At the receiver, the received signals are sampled by the A/D converter and input to the 
synchronization module 16. A correlator is used to detect the beginning of a packet. The 
output of the correlator also comprises the impulse response h t of the channel. The 
correlation is performed between the received input samples and the expected 
synchronization sequence which must be known a priori. 

Note that in one embodiment, the channel estimation assumes a fixed length channel 
denoted by m. If the channel is actually longer, then the resulting estimate will lose some 
information. If the channel is actually shorter, then some processing would be performed 
unnecessarily. Optionally, the channel length may be estimated before the channel estimate 
is performed using signal processing techniques well known in the art. 

In accordance with the invention, the iterative equalizer subtracts the influence of all 
the previous symbols from the current received symbol using hi and the soft decisions 
B n _ v .M n _ m of the previous received bits as expressed in Equation 2 above. Subtracting the 
influence of the previous symbols provides an estimate of the transmitted symbol B n . The 

accuracy of the soft decisions used in the first iteration is improved by performing additional 
iterations. 

In the second iteration, the previously calculated soft decisions are updated so as to 
generate more precise soft decisions which are subtracted from the current received symbol. 
At this point in time, estimates for future symbols (i.e. symbols indexed by n...n+m) have 
been calculated and are subtracted. Note that in the general case, the iterations continue until 
the error difference between the transmitted bits and the received bits converges to a constant 
BER (or less then a predetermined threshold). In the example embodiment presented herein, 
the iterative equalizer performs four iterations. The four iterations are fixed and are 
performed regardless of the error difference. Alternatively, the number of iterations can be 
made variable, depending on any suitable error criteria. 

This method of equalization has the advantage of eliminating the need to perform 
complex calculations to determine the inverse response of the channel This is due to the use 
of the matched filter in determining the soft decisions. 

A flow diagram illustrating the iterative equalization method of the present invention 
is shown in Figure 3. As described previously, the iterative equalization method of the 
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present invention requires knowledge of the impulse response h t of the channel (step 30). 
The channel impulse response is determined using the correlator which functions to correlate 
the transmitted synchronization sequence with the expected sequence. In addition, the 
method uses the previous m-l soft decisions which are stored for this purpose (step 32). 

The influence of all the symbols already estimated apart from a window of 5-1 
neighbors, is subtracted from the current symbol using the channel estimate and the m-l 
previous and m-l future soft decisions (step 34). The soft decision is done using Equations 7 
and 8 in step 35. The error difference is then calculated between the symbols estimated 
during this iteration and those of the previous iteration (step 36). The error difference is 
compared to a predetermined threshold (step 38), and if greater, another iteration is 
performed. The soft decisions for the symbols are updated (step 40) and the method repeats 
steps 34, 36 and 38. The method terminates when the error difference is less than the 
predetermined threshold. Note that alternatively, the number of iterations may be fixed as in 
the example implementation presented hereinbelow. 

A block diagram illustrating the iterative equalizer of the present invention in more 
detail is shown in Figure 4. Note that the implementation described hereinbelow has been 
optimized for compact implementation in VLSI and thus does not explicitly implement the 
general formulas presented hereinabove. Mathematically, however, the implementation 
given below is equivalent to the equations presented above. 

In accordance with the invention, the channel coefficients h t are calculated at the 
beginning of each packet. In addition, they can be updated during the packet if needed. The 
iterative equalizer uses the coefficients calculated for the duration of the packet. The channel 
parameters are calculated using the cross correlation coefficient h 10 27 (Figure 2), 

The function of the equalizer is to remove the channel interference from the received 
signal. Under the second assumption that the time response of the channel for a symbol 
stream is equal to the channel response for each individual symbol, the received symbol is 
reconstructed by subtracting the time response (i.e. impulse response) of all the other symbols 
from the current symbol. 

The iterative equalizer, generally referenced 18, comprises a modified received 
symbol calculator 40, matched filter 42, soft decision calculator 44, quantizer 46, multiplexer 
48, summer 50 and delay 52. The equalizer is operative to generate soft symbol decisions B n 

using as input the received symbol x n samples, channel impulse response h t from the channel 
estimator and h 10 computed from it. 
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First, the modified received symbol is calculated using block 40. The modified 
received symbol calculator 40 functions to subtract the impulse response of previous 
estimated symbols leaving the impulse response of the current symbol y(ri) 41 cleaner (see 
Equation 2). In the second and higher iterations, only a correction on the previous estimate 
needs to be subtracted. The modified received symbol is then correlated by h t using block 42 
(i.e. matched filter operation). The output 43 of the matched filter, denoted by MF(n), along 
with MF(n-l) and CORR(n-2) from the previous symbol (not the previous iteration) and also 
from th e previous iteration (M), is then used to calculate the transmitted symbol B l n _ x 
(where the superscript denotes iteration T). The value MF{n-\) is generated by passing 
MF(n) through a symbol delay line 47 which delays the matched filter output by one symbol 
time (i.e. four clock cycles). Blocks 40 and 42 will now be described in more detail. 

A block diagram illustrating the modified received symbol calculator/matched filter of 
Figure 4 is shown in Figure 5. The modified received symbol calculator/matched filter, 
generally referenced 80, is operative to implement the modified received symbol calculator 
40 and the matched filter 42 blocks of Figure 4. The apparatus 80 comprises a shift register 
array of four rows of sample shift registers 75, each row comprising 69 sample registers 62 
containing 8 bits of resolution each. A timing and control block 74 controls the timing of the 
sample shift register filling process. The arriving input samples 70 are written into input 
registers 68, one input register associated with each sample shift register. 

Note that the length of the equalizer is in this example is equal to 64 samples. The 
ratio between the sample rate and the symbol rate (or bit rate) is in this example is four 
samples per symbol while the number of iterations is also fixed at 4. The equalizer structure 
disclosed herein can be modified by one skilled in the art to accommodate different numbers 
of iterations and different number of samples per symbol, both of them not necessarily being 
equal to each other. 

The samples arrive to the sample shift register and are written therein as follows. The 
input sample vector (i.e. register 68) is filled after the arrival of every four samples 
comprising a symbol. In this example embodiment, each sample comprises 8 bits of 
resolution. At each sample time, switch 70 rotates so as to steer the input sample to the 
appropriate next register. After the input sample vector is filled with four samples, it is 
parallel loaded into the sample shift register via the plurality of switches 72. At each of the 
intermediary samples, each of the four shift registers are shifted to the right. The shift 
registers are shifted in a cyclic manner due to the position of the quad ganged switch 72. 
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During the three samples, the switch 72 is set to cycle the rightmost sample back to the input 
via loop back connections 77. 

The equalizer is operative to output a soft decision of the current symbol every 
symbol time (every four sample clocks). The equalizer is operative to keep each sample 
within the sample shift register for four cycles through its respective shift register row before 
it is overwritten with a new sample. Each shift register row is of length 69 samples and a 
new sample enters the row every four sample clocks. Thus, on the fifth cycle through the 
row, a sample is overwritten with a new one. In general, the length of the row must be one 
more than a multiple of the number of samples per symbol. 

For illustration purposes, the first 16 input samples, labeled Xi through X 16 , are 
shown in their respective sample locations in the sample shift register array. Every four 
sample clocks, the contents of the input sample register are loaded into the sample shift 
register array. The first symbol, for example, comprises sample locations labeled Xi, X 2 , X 3 
and X 4 . The sample shift register array thus holds 64 samples for each iteration at any one 
time. 

As the contents of each row is cyclically shifted right, the channel estimate h is 
multiplied, via multipliers 64, with a CORR value representing the new symbol decision 
minus the decision for the same symbol calculated during the previous iteration. The result 
of the multiplication is subtracted from the input sample via adders 66. The channel estimate 
comprises 64 coefficients. This represents a channel length m of 16 symbols with each 
sample comprising four samples. It is appreciated that the equalizer of the present invention 
is not limited to a channel length of 16 and can be modified by one skilled in the art to 
accommodate other channel lengths as well. 

Each row comprises 16 multipliers and 16 adders. The top row uses coefficients h 63 , 
h 59 ...h 3 ; while the next row uses coefficients h 62 , h 58 ...h2; the following row coefficients h 6h 
h 57 ...hu and the bottom row coefficients h 60 , h 56 ...h 0 . Thus, at each sample time, the 
influence of the previously processed symbols is subtracted out. Sixteen subtractions are 
performed for each row. Each sample, however, is shifted through its row four times 
resulting in four symbol subtractions per input symbol. Each such subtraction is timed such 
that it is carried out on the sample belonging to next iteration. Likewise, the CORR 
calculation is timed such that it is available at the right time in the right iteration. In this 
fashion, the sample shift register array implements the iterative subtraction described supra. 

It is important to point out that the position in each row indicated by reference 
numeral 73 does not have a multiplier/adder associated with it. Once the correction impulse 
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response of the previously processed symbol has been iteratively subtracted, the resulting 
modified received symbol 41 is effectively passed through a matched filter via multipliers 84. 
A total of 64 taps of the sample shift register array contents are connected to multipliers 84 
which are operative to generate the product of the channel response and the modified 
5 received symbol. A tap is connected to every fourth register in each row. Sixteen multipliers 
are associated with each shift register row, wherein coefficients h 63 , h 59 ...h 3 are associated 
with the first row; coefficients h 62 , h 58 ...h2 are associated with the second row; coefficients 
heu h 5 7."hj are associated with the third row; and coefficients h 60 , hs 6 ...h 0 are associated with 
the fourth (bottom) row. The apparatus 80 thus effectively implements Equation 4 on the 

10 current content of the modified received signal. 

The outputs of the multipliers may be stored in product registers (not shown) Summer 
88 is adapted to add all 64 product results so as to yield a matched filter soft output 43 for the 

O current symbol, denoted MF(n). 

m Note that since it is convenient to denote the last computed matched filter output by 

;t5 the index n, instead of n+1, the index n will be shifted by 1 in the following with respect to 
previous section. With reference to Figure 4, the output 43 of the matched filter is input to 
3 the soft symbol calculation block 44. The block utilizes the channel estimate coefficient h 10 , 
™ the MF ( W ) from the matched filter and the previous MF(n-l) and the previous CORR value to 
£ calculate the transmitted symbol B' n _ x . The calculations performed by the block 44 are 
20 expressed mathematically by the following 

m 

7 (9) 
where k is equal to number of samples per symbol, i.e. 4, and 
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At this point it is important to realize that the contents of the sample shift register are 
influenced by the operations in the previous iterations (or steps). From Equation 10 we 
obtain 



B = 



DET 



' i Z n ' ^00 Z n-l ' ^10 ] 
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Where DET is the determinant of the 2x2 matrix in Equation 10. Note that by way of 
experimentation by the inventors, it was found that division by the determinant DET may 
optionally be omitted without sacrificing performance. In addition, it was found that 
5 normalization by dividing by h 00 is required to achieve adequate performance. 

The values z n and z n -i are obtained from the matched filter output MF. It is important 
to note that the values z n and z n -i are calculated in order to compensate for the terms that were 
removed during previous operations. The effects of these terms are calculated and placed 
back into the expressions for z n and z n -u The calculation is as follows 

10 = MF n _ x + ^ . CORR n _ 2 + h m B^ n (12) 

This removes the effect of the last correction, where CORR(n-2) is the correction applied on 
the previous symbol and need to be compensated 

% Zn=MF n +B?hto + K*-B% (13) 

J I This effectively removes the influence of the previously estimated B„. lt B n which was 

|p removed from the sample shift register in the previous iteration and thus from the value MF. 

Note that the inventors not only detected no loss in performance, but actually a gain was 

% f noticed, when the term B n h m was removed from Equation 13. The reason for this is that we 

H want to estimate B n _ x so we would like to use its influence on the received signal. 
T: The value B n _ } is needed, thus we compute 

20 S n-l = -~~ * {Z n * (-V + Z n ^ • /*oo ) (14) 

The result B n ^ is then (optionally) input to the soft decision function 45. Note that some 

optimization in the number of calculations can be obtained by substituting Equation 12 and 
13 into 14. 

The output of the soft decision function is, for implementation proposes, quantized 
25 into one of a finite number of levels via quantizer 46. In the example embodiment presented 
herein, 5 levels are used {-1, -1/2, 0, 1/2, 1}. The quantized and soft limited version of the 
symbol B n _ x is then subtracted from the corresponding symbol from the previous iteration via 
subtractor 50. Thus, the previous estimate of the same symbol is subtracted from the current 
estimate. Ideally, the value of CORR should decrease to a minimum. The difference is 
30 denoted by CORR and is used by block 40 in removing the influence of the correction in the 



12407.0014 



22 



value of B n _ x and is also used by block 44 to calculate B n ^ . Note that for each symbol, four 

values of B n _ x are calculated, each during a different iteration of the equalizer. 

The output of the quantized B l n _ { forms the output of the equalizer and is input to the 

subsequent stage, i.e. the de-interleaver or decoder (if interleaving is not employed). Note 
5 that for each symbol, four values of B l nA are generated. Only the final value, i.e. the fourth 

value when i=4), calculated after the fourth iteration for the symbol, is output to the next 

processing stage. The first, second and third values are ignored. 

In addition, the output of the quantization B l n _ x is delayed one complete iteration delay 

via iteration delay line to yield 52. The iteration delay 52 functions to compensate for 
10 the sample shift register delay encountered in the shift register rows within blocks 40 and 42 
O b y imposing a delay of about 68 clocks. The output of the iteration delay is input to 
:jf calculation block 44 and to one of the two inputs to multiplexer 48. The multiplexer 

functions to pass the Rvalue during the second, third and fourth sample clock cycles and 
l;:; zero during the first clock cycle since no previous estimated symbols are generated until the 
|S second iteration. The output of the multiplexer is subtracted from the quantized version of 
O K-x via subtracter 50 to yield the CORR value. 

r\ Each time a symbol traverses the sample shift register four times, a new correction 

% factor CORR is calculated and a new soft symbol is generated. A symbol returns for 
I- updating and re-decision every 69 clocks, therefore the last B l n ^ used for the current symbol 
20 is used again after an iteration delay of 68 clock cycles to calculate the new B^ value for 
the next cycle. After four iterations in the equalizer, the soft decision for B l n __ x is assumed to 
converge to a constant BER. Note that depending on the application this may not be a valid 
assumption. Therefore, the number of samples per symbol, number of clock cycles per 
symbol and number of iterations may be fixed at a value other than 4, and are not necessarily 
25 equal to each other. 

In an alternative embodiment, the equalizer is constructed using concatenated shift 
registers rather than the cyclic shift register structure of Figure 5. For this embodiment, 
however, additional hardware is required including duplicate multipliers and adders. 

Note also that the output of the equalizer may (1) be input to the input of ECC 
30 decoder if one is used in the receiver, (2) be input to a de-interleaver before the ECC if both 
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are used or (3) comprise the output of the receiver directly if no interleaving or ECC is used. 
Further, the output of the equalizer may be taken before or after the soft decision function 45 
(Figure 4). Due to coarse quantization used in the equalizer example above in order to save 
complexity, the decoder can be modified to use a larger number of bits in the soft input. 

5 Computer Embodiment of the Iterative Equalizer 

In another embodiment, a computer is operative to execute software adapted to 
perform the iterative equalization method of the present invention. A block diagram 
illustrating an example computer system adapted to perform the iterative equalization method 
of the present invention is shown in Figure 6. The system may be incorporated within a 
10 communications device such as a receiver or transceiver, part of which is implemented in 
software. 

O The computer system, generally referenced 140, comprises a processor 142 which 

may be implemented as a microcontroller, microprocessor, microcomputer, ASIC core, 
m central processing unit (CPU) or digital signal processor (DSP). The system further 
fi$ comprises static read only memory (ROM) 144 and dynamic main memory, e.g., random 
° access memory (RAM) 154, all in communication with the processor. The processor is also 

in communication, via bus 146, with a number of peripheral devices that are also included in 
JE the com P u ter system. An A/D converter 156 functions to sample the received signal output 

of the Rx front end 152 coupled to the channel 150 (e.g., the power line, twisted pair, etc.). 
|f Samples generated by the processor are input to the front end circuit via D/A converter 157. 

The front end circuit comprises PHY circuitry, including receiver, transmitter and channel 

coupling circuitry. 

One or more communication lines 158 are connected to the system via I/O interface 
160. A user interface 162 responds to user inputs and provides feedback and other status 

25 information. A host interface 166 connects a host device 164 to the system. The host is 
adapted to configure, control and maintain the operation of the system. The system also 
comprises magnetic storage device 148 for storing application programs and data. The 
system comprises computer readable storage medium which may include any suitable 
memory means including but not limited to magnetic storage, optical storage, semiconductor 

30 volatile or non-volatile memory, biological memory devices, Jaz, Zip, DAT or any other 
memory storage device. 

The iterative equalization software is adapted to reside on a computer readable 
medium, such as a magnetic disk (e.g., Zip, Jaz) within a disk drive unit or magnetic tape 
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means (e.g., DAT). Alternatively, the computer readable medium may comprise a floppy 
disk, Flash memory card, EEROM, EPROM or EEPROM based memory, bubble memory 
storage, ROM storage, etc. The software adapted to perform the iterative equalization 
method of the present invention may also reside, in whole or in part, in the static or dynamic 
5 main memories or in firmware within the processor of the computer system (i.e. within 
microcontroller, microprocessor, microcomputer, DSP, etc. internal memory). 

In alternative embodiments, the method of the present invention is applicable to 
implementations of the invention in integrated circuits, field programmable gate arrays 
(FPGAs), chip sets and especially application specific integrated circuits (ASICs), including 
10 wireless implementations and other communication system products. 

Example C code that implements the iterative equalizer of the present invention is 
presented below in Program Listing 1. The program code may be stored in magnetic 
O memory, disk drive, tape unit, firmware volatile or non-volatile memory for execution by the 
f ^ computer system described above. 

I? 

ft! Program Listing 1 : Iterative Equalizer C Code 

// The assumption here is h00=l (i.e. the channel is normalized) 
- // h 1 0 is precomputed 

O // T is array of soft decision slopes, different slope per iteration. 

2Q II L is the channel length in samples 

O // IT is number of iterations 

: :J // the code assumes one sample per symbol for simplicity, which may be complex. 

' " int Decode(complex x) // x is one sample input to the equalizer 

25 { 

int a, it, i; 

complex Delayed_B, MF, t, B_out; 
complex correct, Z_out, Delayed_out; 
complex soft, old; 

30 

Z[0][0]=x; 
Delayed_B=0; 

for (it=0; it<IT; it++) { // it runs on the iterations 
MF=0; 

35 for (i=0; i<L; i++) MF Z[it][i]*conj(h[L-l-i]); // implement the matched filter 

// last_MF, last_correct are MF, correct of previous symbol respectively 
t=-hlO*s+lastJMF[it] + conj(hlO)*last_correet[it] + Delayed_B*(l-hl0*hl0); 
// Delayed_B is B delayed by a full shift register length, so it is the u of the previous 
40 iteration 

B=f(t, T[it]); // pass through the soft decision function 
correct=Delayed m B -B ; 
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// feedback 

if (I<L) correct=0; // I is the bit number, and if I<L it means the shift register is not 
full, so correct is not valid 

// perform the update of the shift register contents. 
5 Z_out=Z[it][L] + h[0]*correct; 

for (i=L; i>l; i-) Z[it][i]=Z[it][i-l] + h[L-i+l]*correct; 
Z[it][l]=Z[it][0]; 

if (it > 0) Z[it][0]-last_lteration_out; // output of last iteration delay line 
1 0 last_Iteration_out=Z_out; 
// delay of B 

Delayed_B=DLYU[it].inout(B); // DLYU is an array of IT delay units delayed by 
L+l each. 

15 

last_s[it]=s; 
last_correct[it]=correct; 

£3 } 

ip 

return B; 



jr: Further, it is appreciated that one skilled in the computer and electrical arts can 

25 implement the code in Program Listing 1 as a combination of hardware elements in a 
^ hardware implementation without departing from the scope of the invention. In this case, the 
J-J each iteration is implemented using a separate shift register. Note that the number of 
O multipliers in a hardware implementation of the software program listing is increased since 

their operating speed is now lower. Using multipliers as fast as in the iterative hardware 
30 embodiment described above (i.e. single shift register case) permits four multiplies using a 

single multiplier. In this case, however, the same number of multipliers are needed and thus 

the same number of MIPS is needed for both implementations. 

Further, the invention may be implemented in the alternative in a software 

embodiment wherein the data samples are stored in a block buffer and the processing 
35 performed on a block by block basis rather than on a continuous basis as in the iterative 

hardware implemented described above. In addition, the processing need not be performed 

sequentially symbol by symbol but may be performed in any order desired. 

It is intended that the appended claims cover all such features and advantages of the 

invention that fall within the spirit and scope of the present invention. As numerous 
40 modifications and changes will readily occur to those skilled in the art, it is intended that the 

invention not be limited to the limited number of embodiments described herein. 
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Accordingly, il will be appreciated that all suitable variations, modifications and equivalents 
may be resorted to, falling within the spirit and scope of the present invention. 
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